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SUMMARY 

A study was made of the gust response of an annular turbine cascade using a two-dimensional Navier 
Stokes code. The time-marching CFD code, NPARC, was used to calculate the unsteady forces due to the fluid flow 
The computational results were compared with a previously published experimental data for the annular cascade 
reported in the literature. Reduced frequency, Mach number and angle of incidence were varied independently and 
the gust velocity was sinusoidal. For the high inlet velocity case, the cascade was nearly choked. 


NOMENCLATURE 


A,B amplitudes 
C h blade chord length 

M Mach number 

N number of cylinders 

p static pressure 

S blade spacing 

T temperature 

U t axial inlet velocity 

t time 

u»v axial and tangential velocity components 

x,y axial and tangential coordinate directions 

Greek Symbols: 

a phase angle 

(} incidence 

<t> flow coefficient 

p density 

a solidity (C h /S) 

co h angular velocity 

Subscripts: 

1 inlet conditions 

2 conditions downstream of rotor 
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rotor conditions 
radial conditions 
blade properties 


w 


r 

h 


INTRODUCTION 

Periodic wake-passing from upstream blade rows in the gas turbine flow field strongly influences the 
unsteady forces on the surlaces of the downstream blades. Manwaring and Wisler ( 1993) conducted a series of 
experiments involving a gust inlet velocity in a large speed compressor and turbine rig and provided aerodynamic 
forced response measurements on downstream stators. Kurkov and Lucci (1997) provided experimental data on a 
gust response of an annular turbine cascade in a transonic flow regime. They varied the reduced frequency. Mach 
number and incidence angle independently. The mean flow was documented by measuring blade surface pressures 
and cascade exit flow. 

The unsteady aerodynamics generated by a strong blade row coupling and the interactions due to the highly 
unsteady flow induce vibrations, which can lead to high cycle fatigue. In order to predict blade life or design blades 
for longer high cycle fatigue life, accurate predictions of the multiblade row unsteady aerodynamics are crucial. 
Turbomachinery aeroelastic analyses have been carried out at NASA Glenn Research Center for the past decade. 
Kielb and Chiang ( 1992) reviewed research related to the forced response analyses of turbomachinery. Hall and 
Verdon (1990) provided an analysis lor gust response lor cascades. Giles ( 1988) computed the interaction of an 
unsteady wake with a rotor. Reddy and Srivastava (1994) presented an analysis of a rotor- stator stage interaction 
by solving Euler equations. Reddy and Srivastava (1996) extended their previous work to include the effects of 
incoming gust. Barter et al. (2000) studied the interaction effects in a transonic turbine stage by means of a three- 
dimensional. viscous, time-accurate computational code. Their results showed that interaction effects must be taken 
into account in order to accurately predict the unsteady loading on the upstream blade row. Haldeman et al. (2000) 
reported a combined experimental and computational effort quantifying unsteady airfoil interactions. 

The present work was undertaken in order to study the fluid/structure interaction for the turbomachinery 
blade configuration using the state of the art CFD techniques. The time-marching CFD code, NPARC was used to 
calculate the unsteady forces due to the fluid flow . The computational results are compared with the experimental 
data of Kurkov and Lucci ( 1997). 


GOVERNING EQUATIONS 


The differential equations used are the Reynolds mass-averaged Navier-Stokes equations inside a rotating 
blade passage and may be written as: 


p, ; +(pU,)„ = 0 

(P Ui )>, +(p UiVj ),,• + 2p t ijk L1jU k = - Pn + pfe/,, + Ujj - 1 U u 5,j j - p^r ; 


'i +F i 


(P e~p\, +(P Vi e \j = 


(£),/•; V -Kb ♦ v,F, + u u , 


p-pRT 


where U t = mean velocity, u j = fluctuating velocity, e = total energy, £2- = angular velocity. For stationary blades, 
w'e set £2- = 0. It is well known that neither the conventional mixing length type turbulence model nor any standard 
two-equation type turbulence model describes turbulence stresses properly in the region behind a shock wave or in 
the separated flow' regions. Several recent studies have shown that significant improvement can be achieved when 
the standard tw'o-equation model is modified to include the low Reynolds number effects. The following are the 
transport equations for the turbulent shear stress: 
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Standard expressions were used for |i^ t and Cs. These details are omitted in ordei to conserve space. 
The mean flow and turbulence equations were integrated in time using a fully coupled approximately factored 
implicit backward Euler method. 


CODE DESCRIPTION 

T\vo-dimensional time-marching How simulations were performed wath NPARC version 3.1 code. NPARC 
solves the Euler or Navier-Stokes equations in conservation law form on a multi block body fitted grid system. The 
fluid is assumed to be ideal and Newtonian and the Fourier heat conduction law' is assumed for heat transfer. The 
flow can be assumed to be laminar, turbulent or inviscid, A variety of turbulence models, including the k-e model 
can be selected. The flow simulations were performed in two steps. First, a converged steady state flow solution 
was achieved. In the second step, a convective gust perturbation of known strength was introduced at the upstream 
boundary and allowed to propagate down the cascade. The steady state simulation to initialize the flow field was 
done with the approximate factorization algorithm ({SOLVE = I) using local time stepping (IVARDT = 2) with 
DTCAP = 1. Typically, about 2500 time steps w>ere required to reduce the L 2 norm to about 10~ 7 . For the time 
accurate computations, a five step Jameson algorithm second order accurate in time was selected ([SOLVE — 5). The 
time step throughout the grid block was set to a constant equal to the CFL limit specified by DTCAP (IVARDT = 4 
and DTCAP = 0.3) at the location of maximum change in the flow variables. 


GRID GEOMETRY DESCRIPTION AND BOUNDARY CONDITIONS 

Figure 1 shows a single grid block used for computations. The airfoil upper and low'er surfaces are located 
along lines BC and FG, respectively w'here solid wall boundary conditions are prescribed. Lines AE and DH repre- 
sent inflow and outflow boundaries, respectively. The flow incidence angle is specified and the backpressure is 
adjusted until the average Mach number along the inflow boundary AE matches a specified value. Periodicity is 

imposed between lines AB and EF and lines CD and GH. 

Computational grids were generated using the GRIDGEN2D grid generation program. A uniform grid 
spacing w ; as used in each coordinate direction. A grid refinement study was completed to ensure that the computed 
results were independent of the grid density. A grid of 41x1 18 was selected lor the cascade simulation. There were 
60 points on the airfoils, 29 points between the inlet and leading edge and the rest between the trailing edge and 
exit plane. 

For treating the velocity perturbations, namely, the gust for the present problem, the wake profile was fitted 
to experimental velocity profile at the inlet in the following form: 


1 1 + ASin 

N w \ 

( y ^ 

l + Cti \ 

1 


U j 

1 Ml 

l + BSin 

nJ 

f v ) 

+ «T } 

1 

H 1 

U J 

' -Jj 


In the above equations, A and B are the amplitudes; Uj the axial inlet velocity and (Xj and <x 2 the phase angles. 
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DESCRIPTION OF THE EXPERIMENTAL SETUP 


The air entered the belimouth and passed through the annulus to the rotor consisting of radial pins, either 
I 7 or 0.476 cm ( I /8 oi 3/16 in. ) in diameter. The wakes from the pins were eoiivectcd by the main flow' into a 
turbine cascade. The rotor pins and the turbine cascade are shown in figure 2. The rotor induced a swirl of 3 to 5° 
inti) the mean flow. The cascade turning angle was 6 1 °, the blade chord angle relative to the axial direction was 38° 
and the distance between the rotor and the cascade leading edge was 3.9 axial chords. The chord length was 4.75 cm 
and annulus outside and inside diameters were 40.64 and 27. 1 3 cm, respectively. The cascade consisted of 23 blades 
that extended in the radial direction spanning the annulus. The power was delivered to the rotor through a long shaft 
that extended vertically through the inner pipe of the annulus. The shaft exits the outer pipe at the 90°-bend close to 
the floor level and is terminated by a belt-driven pulley, which transmits the power from the electric motor. 

The number of pins on the rotor was varied among 24, 12 and 6. The pin diameter for the 24-pin configura- 
tion was 0.3 17 cm and for 12 and 6-pin configurations, it was 0.476 cm. The pin diameter was increased for the last 
two configurations in order to reduce the number of harmonics in the wake. The reduced frequency 

(0 = ■ to r < 'h _ n 77 ^ ^ » 

{/, ' 0 was varied by changing the number of pins on the rotor, N w . Here, to r is the shaft angular 

velocity. c h the blade chord, U| the axial inlet velocity and 6 the flow coefficient. Subscript 1 denotes the station 
downstream oi the belimouth and upstream of the rotor. The reduced frequency values were 10 , 5 and 2.5 for the 
available values of N w . The inlet Mach number was either 0.2 or 0.27. At the higher inlet Mach number the flow in 
the blade passage throat was near choking at the midspan radius. The flow coefficient was kept constant at 0.9 and 
the shaft speeds for the two Mach numbers were 4290 and 5790 rpm. The Reynolds numbers corresponding to low' 
and high inlet velocities were 2.7xl0 5 and 3.6x10 s . respectively based on the blade chord. Positive and negative 
values of incidence correspond to two directions of rotation. 

The inlet condition was defined by the barometric pressure, inlet temperature at the belimouth inlet and 
several static pressure measurements on the inside and outside wall of the annulus downstream of the belimouth. 

The cascade exit flow' w'as surveyed at the blade midspan radius over a circumference in excess of one blade pitch 
with a five-hole probe. This station was one axial chord length downstream of the cascade trailing edge. The loca- 
tions of the steady state and dynamic ports on the instrumented blades are shown in figure 3. Kurkov and Lucci 
(1997) have documented in more detail the measurements taken. 


RESULTS AND DISCUSSION 

Figure 4 displays the steady state results for the pressure coefficient distribution for inlet Mach number 
M| = 0.2 and inlet angle of attack of 4.59°. Figure 5 shows the steady state results for the pressure coefficient distri- 
bution tor inlet Mach number M | = 0.27 and inlet angle of attack of 5. 17°. The computed results agree with the 
experimental data quite well in both cases. The spikes in the pressure coefficient distribution in the vicinity of the 
trailing edge may be due to the definition of the airfoil surface coordinates in generating the grid. 

Figures 6 to 9 show results for the distribution of the amplitude of the unsteady pressure coefficient (per- 
centage of the inlet dynamic head) for inlet Mach number M, = 0.2 and reduced frequency of 10 based on chord 
length. Figures 10 to 13 show results for the distribution of the amplitude of the unsteady pressure coefficient for 
inlet Mach number M , = 0.27 at the same reduced frequency. The computed results agree with the experimental 
results within 20 percent error over most of the blade length, both on the suction and the pressure surfaces. How- 
ever, in the vicinity of the leading edge of the blade, the discrepancy seems to be higher. 

Figures 14 to 17 show results for the distribution of the phase angle for inlet Mach number M, = 0.2 and 
reduced frequency of 10. Figures 18 to 21 show the same results for inlet Mach number M, = 0.27. The quantitative 
agreement between the computed and measured values of the phase angle is poor although qualitatively they are 
showing similar trends. 
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CONCLUDING REMARKS 


In this paper, we have considered an annular cascade configuration subjected to unsteady gust inflow 
conditions. The gust response calculation has been implemented into the time marching CFD code. NPARC. The 
computed steady state results tor the pressure distribution demonstrated good agreement with experimental data. 
The computed results for the amplitudes of the unsteady pressure over the blade surfaces agreed with experimental 
data within 20 percent error margin, except in the vicinity of leading edge. There exists considerable quantitative 
discrepancy between the computed and experiments results for phase angle, although the trends arc similar. 
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Figure 2. — Experimental setup. 
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Figure 4.— -Steady state pressure prediction (Mi = 0.2, i « 4.59°). 



Fraction of chord 

Figure 5. — Steady state pressure prediction (M-j = 0.27, i * 5.17°). 
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Ntude 



Figure 6— Amplitude of the unsteady pressure coefficient distribution 
on suction surface (Mi = 0.2, i = 3.24°). 



Figure 7.— Amplitude of the unsteady pressure coefficient distribution 
on pressure surface (Mi = 0.2, i = 3.24°). 
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Figure 8.— Amplitude of the unsteady pressure coefficient distribution 
on suction surface (Mi = 0.2, i = 3.24°). 



Distance along chord 

Figure 9. — Amplitude of the unsteady pressure coefficient distribution 
on pressure surface (Mi = 0.2, i = 3.24°). 
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Figure 1 0.— Amplitude of the unsteady pressure coefficient distribution 
on suction surface (Mi = 0.27, i = 2.98°). 



Figure 1 1 .—Amplitude of the unsteady pressure coefficient distribution 
on pressure surface (Mi = 0.27, i = 2.98°). 
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Distance along chord 


Figure 1 2.— -Amplitude of the unsteady pressure coefficient distribution 
on suction surface (Mi = 0.27, i = 2.98°). 



Distance along chord 


Figure 1 3. — Amplitude of the unsteady pressure coefficient distribution 
on pressure surface (Mi « 0.27, i = 2.98°). 
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Phase angle 3! Phase angle 



jure 1 4 . — phase angle distribution on suction surface (Mi = 0.2, i = 3.24°). 



Figure 1 5.— Phase angle distribution on pressure surface (M t = 0.2, i = 3.24°) 
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Distance along chord 

Figure 16. — Phase angle distribution on suction surface (M^ = 0.2, i = 3.24°). 



Distance along chord 


Figure 17. — Phase angle distribution on pressure surface (M-j = 0.2, i = 3.24°). 
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Phase angle g Phase angle 



ure 18, — Phase angle distribution on suction surface (M^ = 0.27, i - 2.98°) 



Figure 19. — Phase angle distribution on pressure surface (M-j = 0.27, i = 2.91 
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